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Abstract 

We present a quantum nonlinear model of two-wave mixing in a lossless photorefractive 
medium. A set of equations describing the quantum nonlinear coupling for the field operators 
is obtained. It is found that, to the second power term, the commutation relationship 
is maintained. The expectation values for the photon number concur with those of the 
classical electromagnetic theory when the initial intensities of the two beams are strong. 
We also calculate the quantum fluctuations of the two beams initially in the coherent state. 
With an appropriate choice of phase, quadrature squzzing or number state squeezing can be 
produced. 

PACS numbers:42.654.Hw, 42.50.Dv, 42.50.Lc 


1 Introduction 

The photorefractive effect in electro-optic crystals, a phenomenon in which the local index of 
refraction is changed by the spatial variation of light intensity [l], has been studied extensively for 
its potential in many applications. The fundamental process may be described as follows. When 
the crystal is illuminated with a spatially modulated intensity pattern, free carriers (for example, 
electrons) are nonunifromly generated due to the photoionization of impurities (generally, which 
may be doped). The impurities that can be ionized and provide free carriers are called donors. 
Once these donors are ionized they can serve as trap sites which capture electrons. The electrons 
can be transported by diffusion or drift and become trapped at these sites. The trapped electrons 
can then be re-excited except for those in the dark region. Thus a space-charge separation is 
created, which leads to a space-charge field. Such a field induces a change in index of refraction 
via the Pockels effect (linear electro-optic effect), creating an index grating. The presence of 
such an index grating will in turen affect the propagation of these beams. Crystals such as 
LiNb0 3 , BaTi0 3 , SBN, BSO, GaAs and InP, are efficient media for the generation of the 
photorefractive effect with relatively low intensity level (eg., lW/cm 2 ). 

Many different nonlinear optical phenomena in photorefractive media have been studied. 
These include wave mixing, phase conjugation, self-oscillation, photorefractive resonance, etc. 
The fundamental photorefractive process is two-wave mixing (TWM), in which two beams in- 
tersect inside a photorefractive midia. A stationary index grating is formed which is sptaially 
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shifted 7r/2 relative to the intensity pattern. Such a spatial phase shift leads to nonrecipro- 
cal energy transfer when these two beams propagate through the medium. The basic classical 
electromagnetic theory explaining the nonlinear interaction involved is already well established. 
Much attention has been focused on its applications including photorefractive resonatore, non- 
reciprocal transmission windows, self-pumped phase conjugators, laser beam clean-up, optical 
interconnection, etc. Although a number of cases of TWM have been analysed, a quantum 
theory is not available and photorefractive non-classical effects have not been discussed. In this 
paper we present, to our knowledge for the first time, a quantum treatment of two-wave coupling 
in a lossless photorefractive medium. 

2 Quantum model of photorefractive TWM 

A typical geometry for studying two-wave mixing is shown in Fig.l. Under certain circum- 
stances, two beams of light can interact in a photorefractive crystal in such a manner that 
energy is transferred from one beam to the other. This process is also known as two-beam 
coupling. The signal and pump waves, of amplitudes A t and A, respectively, interfere to form 
a nonuniform intensity distribution within the crystal. Due to the nonlinear response of the 
crystal, this nonuniform intensity distribution produces a refractive index grating within the 
material. However, this grating is displaced from the intensity distribution in the direction of 
the positive (or effective electrooptic coefficient) crystalline c axis. As a result of this phase 
shift, the light scattered from A p into A, interferes constructively with A s , whereas the light 
scattered from A , into A p interferes destructively with A p , and consequently the signal wave is 
amplified whereas the pump wave is attenuated. 


attanuated Fig.l. Typical geometry for study- 

ing two-beam coupling in a photore- 
fractive crystal. 

amplified 

► z 

An ideal quantum model for degenerate two- wave mixing may be constructed as follows. 
Consider the effective interaction Hamiltonian 

H tff = hx'A^{A^B + B^A)B + h.c. (1) 

where x' * s the effective interaction coefficient for the nonlinear process, A and B are the Boson 
operators for two modes with frequency w 0 = u> h . Factor (A^B+B^ A) represents the interference 
of two modes [2]. The TWM can be understood from the following physical picture. Mode A is 
generated accompanied by the annihilation of mode B, due to scattering from the grating induced 
by the interference. In other words, mode B is “scattered” by the grating in the direction of beam 
A, to yield mode A, which is responsible for the energy coupling. The Heisengerg equations 
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of motion for the field operators A and B may be easilly obtained from H eff . Making the 
conversion z = vt for propagation along the z-axis at a velocity v, we can write the equations 
as 

— = -2i X A^B 2 -t{ X + x)AB^B (2) 

— = -2ix' A 2 - i{x + x*)BA^ A ( 3 ) 

where x = x'/ v ■ We find the field operators satisfy the Boson commutation rules. From the 
equations of motion for the photon number operatore N a and TV* 


dNa 

dz 


= -2i X A^ 2 B 2 + 2ix* A 2 B^ 2 


(4) 


^ = 2i x A^ 1 B 1 - 2.x'A‘B^ (5) 

we can show that the total photon number is constant throught the process. In the short path 
approximation, the solutions of Eqs.(2) and (3) for the field operators with expansion up to the 
quadratic (x 2 ) 2 term is 

A(z) = a — 2i{xz)cfib 2 - iz{x + x'*) 0 ^ 

- I Xz I 2 [(4at a 2 + a)bh - ( 2 at - a)tf 2 b 2 + ( a t 2 b 2 + 6 t 2 a 2 )a + 2 at a 2 ] ( 6 ) 
-z 2 {x 2 {ai 2 ab 2 + b^ 2 b 2 a/2 + 6 + 60 / 2 ) + x* 2 (a 3 6+ 2 + a 6 + 2 6 2 / 2 + abh/2)} 

B(z) — b - 2 i(x*z)b^a 2 - iz(x + X*)bala 

- I X2 | 2 [{4bh 2 + 6 )ata - (26+ - 6 )at 2 a 2 + ( 6 + 2 a 2 + a+ 2 6 2 )6 + 26+6 2 ] (7) 

-z 2 {x* 2 {b^ 2 ba 2 + a+ 2 a 2 6/2 + a+a 6 / 2 ) + x 2 ( 6 3 a + 2 + 6 a+ 2 a 2 /2 + 6 a+a/ 2 )j 

where a and b are the input field operators, respectively. It may be seen that, to the quadratic 
term, the field operators still satisfy the commutation relation 


{A(z),A^{z)] = l, [B{z),B^(z)\ = 1 ( 8 ) 

In order to test our quantum model for photorefractive TWM, we may derive the expectation 
values for the photon number in each beam and verify if the result of the quantum calculation 
concur with those of the classical electromagnetic theory. When the two beams are initially in 
the coherent state | a) and | ft) with 


we obtain 


a =| a 

| exp(iS a /2) 

(9) 

0=101 

exp(i6i,/2) 

(10) 

x=l x 1 

exp(i<j>) 

(11) 

| a | J +4 sin{<j> + S b - S a ) 

+8 | X 2 pi <*0 1 ! u 0 1 1 - 

Ix*l|a 0 l ! +4 | X 2 I 1 (I0| 4 -M<) 
| a | 2 )[1 + c os<ftcos(4> + 6t - 8z) \ 

(12) 
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(N b )= \a\ 2 + \P \ 2 -(N a ) (13) 

where <j>, S b and 6 a are phase angles, depending on the initial condition. We take 6 b — S a = tt 
and write I a =| a | 2 , I b =| 0 | 2 , then obtain 

(N a ) = I a - 4 sin<j> | I hh + 4 j x^ | 2 (V ~ h 2 ) (14) 


(N b ) = I a + I b — (N a ) 


(15) 


According to the classical electromagnetic theory, the coupled equations for photorefractive 
TWM can be written as [3] 

Ah (16) 


dz 
dl _ 2 
dz 


h + h 
hh 


1 -^- (17) 

h + h V ' 

where I x and / 2 are the intensities of beams 1 and 2, respectively, and 7 is the coupling constant 
with 

27m 


zi mi , . 

7 = - -sin<t> (18) 

Xcosu 

Here n\ is the depth of index modulation related to the electro-optic coefficient, 26 is the angle 
between the two beams inside the medium and <f> is the phase that indicates the degree to which 
the index grating is shifted spatially with respect to the light interference pattern. 

By examining the coupled equations, we note that I 2 is an increasing function of z , provided 
7 is positive. This indicates that the energy is flowing from beam 1 to beam 2. The direction 
of energy flow is determined by the sign of 7, which depends on the orientation of the crystal 
axis. The solution for the intensities h(z) and / 2 (z) are 


h(z) = /i(0) 


h{ 0) + / 2 (0)eip(7z) 


h{z) = /j(0)eip(7z) 


Jo 

Ji(0) + hi® )exp(qfz) 
Jo 


(19) 


( 20 ) 


where /i(0) and / 2 ( 0) are the input intensities of beam 1 and beam 2, respectively, and J 0 is the 
sum intensity with I 0 = h{ 0 ) + / 2 (0). In the short path approximation, the solutions can be 
expanded to the quadratic (72) 2 term as 


h{z) = Ji(0) - 72 


Ji(0)J2(0) 

Jo 


+ h z ) 


; h (Q)7 2 (0) 

2/q 2 


[J|(0)~ Ji(0)] 


( 21 ) 


h{z) = Io-h{z) (22) 

Comparing Eqs ; (l4) and (15) with Eqs.(2l) and (22), we find that the results of the quan- 
tum theory are consistent with those of the classical theory, so long as 7i(0) = I a , / 2 ( 0) = 
/ i) 7 = 4 | x | studio and I a >> 1, I b >> 1. When the input intensities of the two 
beams are strong, the effective Hamiltonian H e fj can give and accurate description of the en- 
ergy exchange phenomenon in photorefractive two-wave mixing, as shown in Fig. 2. We can 
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thus conclude that our quantum model for photorefractive TWM is reasonable and successful. 
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Fig. 2. The intensities of two beams versus the effective interaction length (7 z). 
Dashed curve: the intensities of the classical electromagnetic theory I\ and /2, from 
Eqs.(19) and (20). Solid curve: the quantum average photon number N a and N\>. 
The initial intensities I\ = 10 4 , I 2 = 10 6 , respectively. 


3 Quantum statistic of photorefractive TWM 


To discuss the photon number fluctuations of the quantized field we consider the variance 
(A Nj 2 (z)) or the Fano factor 




(aa7M> 

(AN,(z)) 


(23) 


where (A N 2 (z)) = { N 2 (z )) — (Nj(z)) 2 and j = a,b. To obtain the above expression we need 
to find the expectation values for (AiV a 2 (z)} and (AN b 2 (z)). When the input fields are in the 
coherent state | ct) and | /?), after some tedious calculation we may obtain the expectation values 


(AiV a 2 (z)) = I a -8sin4> \ xz \ I a h + 8 \ xz | 2 I a [2I a 2 I b cos(2(f)) +2I a I b 2 (l + sin 2 <f>) -I 2 -I 2 ] (24) 

(A N b 2 {z)) = I b + 8sin<f> \ X z | I a I b + 8 | X 2 | 2 I b {2l b 2 I b cos(2<t>) +2l b I a 2 {l + sin 2 <f>) - I a 2 - I b 2 } (25) 

where we have taken S b — S a — w. Here <f> — ±| corresponds to the maximum energy coupling 
between the two beams. Eqs.(14) and (15) show that the energy flows from beam A to beam 
B when <$> E [0, 7r], This indicates that A is the pump beam and B the signal beam. Let <f> = 
we rewrite Eqs.(24) and (25) as 

(AiV a 2 (z)) = I a - 8 | xz I Ub + 8 | x* | 2 I a \-2l 2 I b + 4/ a / 6 2 + I b 2 - 7 a 2 ] (26) 

(AAT 6 2 (z)) = A + 8 | x* | /./» + 8 | x* I 2 h\~2I b 2 I a + 4 I b I a 2 + I a 2 - I b 2 } (27) 

The Fano factor plotted against the effective interaction length 7 z is shown in Fig. 3, where 

7 = 4 | x I sin<j>I 0 . We see that the pump mode A can be in a squeezed number state, whereas 

the signal mode B becomes super-Poissonion at the some time. 
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Fig. 3 . Fano factors of beams A and B versus 
effective interaction length 7 z for coherent state 

inputs. Dashed curve: F a , solid curve: F b . The 

initial intensity ratio m = = jT,. F a shows 

the sub-Poissoinan statistical feature. 
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Moreover, the signal mode can never becomes squeezed(at least for our solution expanded to 
the second order). The degree of squeezing in the pump mode depends on the initial intensity 
ratio m (here we define m = If m is large (for example, 100), then the degree of squeezing 

will be very small (in the short path approximation). This is reasonable as the energy coupling 


has little effect on the intensity of the pump modes, so the quantum fluctuations will 
reduced greatly. 

The quadrature phase amplitudes of the two beams are defined as 

not be 


Y — A(*)+At{*) y — 

2 a 2t 

(28) 


X _ gl^+gliil Y, - 

^ 0 2 2i 

(29) 

When the input field are in the coherent state, the field variances may be determined explicitly 
to be 

(A X a 2 (z)) = 

\ + | X z | Ibsin(4> + 6 b ) + 2 \ x* | 2 h\h + I a cos 2 (f) - I a cos6 a ( 1 + cos 2 <f>) 
— (2 I a + I b )cos4>cos(4> + <5&)] 

(30) 

(A Y*{z)) = 

| xz | I b sin((j) + 6 b ) + 2\ xz | 2 I b {I b + I a cos 2 <j) + I a cos6 a (l + cos 2 <£) 
-f (2 I a + I b )cos4>cos((f> + <5(,)] 

(31) 

(A X b 2 (z)) = 

i+ | xz | I a sin(6 a - <f>) + 2\xz\ 2 I a [I a + hcos 2 <f> - I b cos6 b (l + cos 2 <j>) 
— (2 /(, + I a )cos<j>cos(6 a - </>)] 

(32) 

(A Y b \z)) = 

| Xz | I a sin(6 a - <t>) + 2 \ xz \ 2 I a [I a + I b cos 2 <f> -|- I b cos6 b ( 1 + cos 2 <f>) 
+ (2 I b + I a )cos<t>cos(6 a - </>)] 

(33) 

With an appropriate choice of phase, both modes can produce quadrature squeezing. For exam- 
ple, when<£ = nj 2, S b = ir and 6 a = 0 , it is obvious that (AX a 2 (z)) and (AX; , 2 (z)) may be less 


than ^ in the short path approximation. The variances plotted against the effective interaction 
7 z are shown in Fig. 4 . We see that both modes can be in the squeezed state. Furthermore, there 
is strong dependence on which of the input modes is strong. The degree of squeezing in the pump 
mode is great when the initial intensity ratio m is small, as shown in Fig. 4. In the reverse case, 
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if the pump mode A is strong, then the degree of squeezing in the signal mode B will be great. 


(AI a 2 ),(AF a 2 ) 


(AX, 2 ),{An 2 ) 



(a) i b ) 

Fig.4. The variances (AX a 2 ), (A Y a 2 ), {AX 2 ) % and { AY h 2 ) when both input field 
are initially in the coherent state, with m = Aq. Dashed curves: the variances of 
X component, which show quadrature squeezing. Solid curves: the variances of Y 
component. 


4 Conclusion 


In conclusion, we have presented a quantum model of photorefractive TWM, which can well 
describe the energy exchange phenomenon in TWM. A set of coupled mode equations is obtained 
and solved in the short path approximation We have also calculated the quantum fluctuations of 
the two modes and find that when both modes are initially in the coherent state, the pump beam 
can become sup-Poissonian, due to the photon flux in the energy transfer. The same qualitative 
result was also obtained in our previous approach from a set of simplified field equations [4]. 
With an appropriate choice of phase relationship, quadrature squeezing can be produced. 


References 

[1] P. Yeh, Introduction to Photorefractive Nonlinear Optics (Wiley, New York, 1993). 

[2] J. Perina, Quantum Statistic of Linear and Nonlinear Optical Phonomena (Kluwer Aca- 
demic publishers, 1991). 

[3] P. Yeh, IEEE J. Quantum Electron, 3, 484 (1989). 

[4] Tang Zheng, Lii Yajun, Wu Meijuan, Li Shiqun and Wu Ling-An, 6th National Conf. on 
Quantum Optics, Jiangxi, China, Oct, 1994 (in Chinese). 


463 




